Shifts of the nuclear resonance in the vortex lattice in YBa 2 Cu 3 07 
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The NMR and NQR spectra of 63 Cu in the CUO2 plane of YBa2Cu3C>7 in the superconducting 
state are discussed in terms of the phenomenological theory of Ginzburg-Landau type extended to 
lower temperatures. We show that the observed spectra, Kumagai et al., PRB 63, 144502 (2001), 
can be explained by a standard theory of the Bernoulli potential with the charge transfer between 
Cu02 planes and CuO chains assumed. 



I. INTRODUCTION 

Already in 1969 Rickayzen has derived the thermody- 
namic correction to the electrostatic Bernoulli potential 
in superconductors and has proposed -to use this correc- 
tion to study the pairing interactional Capacitive mea- 
surements of the Bernoulli potentiaLhave not shown any 
signal of this correction, however.B~til It turned out that 
the thermodynamic correction is cancelled by the sur- 
face dipole,L] therefore, it can be accessed only by mea- 
surements probing the bulk of the superconductor. For 
a long time, there was no experimental method to ful- 
fill this task so that the thermodynamic correction to 
the electrostatic potential and the related charge trans- 
fer have remained a theoretical topidJTij with no exper- 
imental justification. 

The experiment aimed at the charge transfer in the 
bulk of the superconductor has been performed only re- 
cently. Kumagai, Nozaki and MatsudaEj (KNM) re- 
ported very precise measurements of the nuclear mag- 
netic resonance (NMR) and the nuclear quadrupole reso- 
nance (NQR) spectra for two high-T c materials, a slightly 
overdoped YBa 2 Cu307 and an underdoped YBa2Cu40s- 
Comparing NMR and NQR frequencies they found that 
the magnetic field, present in NMR while absent in NQR, 
causes changes of the electric field gradient (EFG) in the 
vicinity of the nucleus. In this way they first observed 
the charge transfer in the Abrikosov vortex lattice. They 
also pointed out that the conventional theory of the vor- 
tex charge completely fails to explain the experimental 
data giving a wrong sign and too small amplitude of the 
magnetic field effect on the EFG. 

In this paper we discuss the KNM experiment in terms 
of a phenomenological theory of Ginzburg-Landau (GL) 
type proposed in Refs. pd]p2| This theory is extended 
by features necessary to describe a layered structure of 
YBCO, in particular the CuO chains which serve as 
charge reservoirs. Unlike in a conventional superconduc- 
tor, where the charge is transfered on a long distance 
from the vortex core to the region between the cores, in 



the YBCO the dominant role plays a local charge transfer 
between the Cu02 planes and the CuO chains. This local 
transfer will be shown to explain the observed sign and 
amplitude of the magnetic field effect on the EFG. At 
low temperatures, there are two important corrections, 
the space variation of the magnetic field and the small 
but finite disorientation of the axis c in individual grains 
of the measured sample. 



A. Electric field gradient effect in the vortex lattice 

Let us brieflji-recall the experiment of Kumagai, Nozaki 
and Matsuda.113 The 63 Cu atom in the Cu02 plane sub- 
jected to the magnetic field, B || c, with a number of 
holes, TV, has NMR frequencies v\ = jcuB and 



"2,3 = ICuB T (C + AN) . 



(1) 



The first term describes the coupling of the nuclear mag- 
netic momentum to the magnetic field, B = 9.4 T in 
Ref. The nucleus 63 Cu has spin 3/2 and its Zeeman 
coefficient is 7 Cu = 11.312 MHz/T. 

The second term of ((!]) is the EFG effect. The nucleus 
63 Cu has a cigar shape elongated in the direction of its 
magnetic momentum and this electric quadrupole feels 
gradients of the electrostatic field due to a non-spheric 
charge of electrons and ions in its vicinity. The largest 
gradient of the electric field (the principal axis of the 
EFG tensor) is parallel to the axis c. The density de- 
pendent contribution to the EFG is attributed to the lo- 
cal interaction between nucleus with (i-states on the Cu 
site. The experimental fitoll3 to the EFG effect yields 
C = 25 MHz and A = 34 MHz per hole. 

The NQR frequency includes the asymmetry correc- 
tion, fNQR == fox/l + 7 7 2 /3, where r] is the difference 
between gradients in the a and b directions scaled to 
the gradient in the c direction. For the Cu atom in 
the plane of YBa 2 Cu3 07 the a-b asymmetry is rather 
weak, a/I + ry 2 /3 w 1 + 2.6 x 10~ 4 . The frequency 
vq corresponds to the NMR in the zero magnetic field, 
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v o = C + ANq, where No is the number of holes in the 
superconducting non-magnetic state. 

A combination of the resonance frequencies, 



Av = u - - (v 3 - v 2 ) ■ 



(2) 



would be zero if the magnetic field has no effect on the 
number of holes on the Cu sites. Actually, this is the 
case for temperatures above T c . Below the critical tem- 
perature, the shift Av is non-zero indicating changes of 
the quadrupole coupling due to the effect of the magnetic 
field on the hole density. Briefly, Av serves as a sensor 
of the charge transfer in the Abrikosov vortex lattice.113 



B. Charge transfer between chains and planes 

In this paper we evaluate the shift Av. To motivate 
our approach we want to recall the way in which the 
experimental data are obtained. The measurement yields 
distributions F2$(v) of the NMR frequencies, see Fig. 3 
in Ref. [l3|. The frequencies ^3,2, used by KNM in formula 
(Q) are positions of the maxima of ^2,3 iy). Since NMR 
lines have a finite width T, each site of a local frequency 
^2,3 (r) contributes to the spectrum by the Lorentzian 



r 



(v-v^ + r 2 



dr 



( V - Vi (r))* + T*- 



(3) 



The integral runs over the 2D volume, f2 = J dr, covering 
the elementary cell of the Abrikosov vortex lattice. The 
magnetic field B and the density TV are functions of r, the 
frequencies ^2,3 depend on r via (|l|). The experimental 
magnetic field is the mean value, B = (B). 

The finite line width is crucial for the interpretation 
of the KNM experiment. The observed contributions of 
the quadrupole shifts are less than 25 kHz which is much 
smaller than the width of the spectral line V ss 140 kHz. 
In this limit, the maximum of the spectral line is given 
by the mean value and the lowest statistical variations, 



2((v l -(v l )) 3 ) 
T*-6((vi-(vi))*y 



(4) 



The variations are important at low temperatures due 
to a strong compression of the magnetic field in vortex 
cores resulting in a strong space-dependence of the Zee- 
man effect. Close to the critical temperature the varia- 
tions vanish and the maxima of the NMR lines approach 
the mean frequencies, ^2,3 ~ (^2,3)- Close to T c one finds 
from (0) and (|) that 



Av « A(N - (N)) . 



(5) 



According to (||), the shift Av reflects how many holes 
have been removed from the Cu02 plane. The in-plane 
charge transfer contributes only via variations giving, as 
shown below, a negligible correction at all temperatures. 



C. Plan of the paper 

The physical picture we assume here is based on prop- 
erties of the non-magnetic state of YBCO. The super- 
conducting state appears mainly in CuC>2 planes, in 
chains the gap is induced by the proximity effect. The 
transition to the superconducting state shifts the chem- 
ical potential in. the, same manner as in conventional 
superconductors.EH This shift is large in planes, where 
the pairing mechanism is located, while it is negligible 
in chains. In equilibrium the electrochemical potential is 
constant so the unequal shift of the chemical potential 
results in the charge transfer between plains and chains, 
until the equilibrium condition is established by the elec- 
trostatic potential. 

The magnetic field partially suppresses the supercon- 
ducting component what leads to a partial backward 
charge transfer. In the numerical treatment we assume 
both mechanisms, the transfer between chains and planes 
and the in-plane transfer. We show that this backward 
transfer between chains and planes brings the dominant 
contribution to the shift Av. 

In the next section we introduce the Lawrence-Doniach 
model of the layered superconductor. In Sec. Ill we de- 
rive a set of equations for the GL function, the vector 
potential, and the electrostatic potential. The charge 
density is obtained from the Poisson equation adapted 
to the layered structure. In the Sec. IV we evaluate the 
distributions ^2,3 (v) from the frequencies ^2,3(1") given 
by ([!]) with the theoretical values of B-b:) and N(r) ob- 
tained from the extended GL theorylL3. In Sec. V we 
present numerical results and discuss the individual con- 
tributions to the experimental data. We show that the 
approximation (|J) applies close to the critical tempera- 
ture, but at low temperatures the space variation of the 
magnetic field leads to appreciable corrections. In Sub- 
sec. C we compare our results with the model for con- 
ventional (not layered) superconductors and explain why 
the conventional model results in the reversed sign and 
too small amplitude of the line shift. In Subsec. D we in- 
clude a correction for samples made of c-oriented grains 
and show that the deviations from the ideal orientation 
are important at low temperatures. Sec. VI contains a 
summary. 



II. LAWRENCE-DONIACH MODEL 

The model of layered superconductors has been pro- 
posed by Lawrence and Doniachta, see also Ref. [l6[ We 
associate the direction z with the c axis of the YBCO 
and introduce a layer index j and an additional index ± 
which specifies the CUO2 layer position, = jDzL^D pp , 
where D = 11.65 A is the period along the axis c and 
D pp = 3.17 A is the distance between two neighboring 
planes. Positions of layers of chains are z| = jD + \D. 
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In the spirit of the tight-binding description, the vari- 
able z is replaced by the layer index. The wave func- 
tion becomes a function of the layer index, ij)(x,y, z) — > 
if>^' c (r), and the 2D coordinate r = (x, y). 

A. Kinetic energy 



In the chains we treat the condensation energy as negli- 
gible. The free energy thus includes only the contribution 
due to the reduced entropy, = — ^j c T 2 y^l — w -. 

The meaning of symbols is analogous to the case of 
planes. Again, the sum over all layers reads f con = 



The kinetic energy is composed of parts parallel to the 
planes and the Josephson coupling along the c axis. The 
kinetic part of the free energy for the motion in the layer 
is treated within the isotropic effective mass, 



P L = 

J kin 



1 



2m* 



(-iftV-e*A5)^| 



(0) 



where is the in-plane component of the 3D vector po- 
tential A(x,y, 
cal, mX = ml. 



Zj ) . The masses in both planes are identi- 
For simplicity we assume the 2D effective 
mass to* in the chain layer to be isotropic, too. 

The motion between layers is represented by the 
Josephson coupling, 



fi 



'3.+- 
kin 



J D 



kin 



J, 



cp 



1>i 



/kin — ^cp ll>j 



3-1 



(7) 



The phase factors 9 of the off-diagonal terms are given by 
the c-component of the 3D vector potential. For B || c we 
can use a gauge in which the c component of the vector 
potential is zero. In this gauge, all three phases are zero, 
9j l = 0. The total kinetic free energy is a sum over all 

layers, F Un = J drj^j, + /£" )• 



B. Condensation energy 



Following Bardeen0'0 we use the Gorter-Casimir free 
energy to describe the condensation. We assume that 
the energy released by the formation of Cooper pairs is 
located in Cu02 planes. In terms of the wave fupction 
the contribution of planes to the free energy readsll3 



3± 
con 



l 



(8) 



where 7 p is the plane contribution to Sommerfeld's 7 (lin- 
ear coefficient of the electronic specific heat), and 



21^1 



21^1 



(9) 



is the fraction of superconducting holes in the given 
plane. The denominator of @ is the total density of 



pairable holes, n L , = 2\il) L A 



which is a sum of the 



superconducting density and the density of normal holes. 



C. Coulomb interaction 

The total density of holes is linked to the 2D density 
of charge, p L j = en' +p,„ tt , which creates the electrostatic 
potential, 



" ^/|r'-r| 2 + (zj;-zp 2 ' 



(10) 



The Coulomb energy is Tc = | / dr^2- L Pjfj- 

It is advantageous to treat the electrostatic potential 
in the momentum representation, 



p t (k)= J dr^(r)e" lkr 



(11) 



We have moved the superscript 1 into subscript and lifted 
the index j since all equivalent layers are at identical 
potentials, iptj = (p L m . The electrostatic potential then 
reads 



with elements of the inverse Laplace operator 



(12) 



L 



L: 



D e 



D 1 



-kD 



2k 1 



-kD 



-kD„ 



. e -k{D-D p 



L 



2k 1 - c- kD 

jj e -kD cp c -k{D cp +D pp ) 



±c 



2k 



1 



-kD 



(13) 



Here D cp = |(-D — D pp ) denotes the distance between 
adjacent chains and planes. 

When inverted, the relation ( |l2| ) yields the density of 
charge in terms of the electrostatic potential, 



3 t (k) = e£>^L tt /^/(k). 



(14) 



In our approach the electrostatic potential is evaluated 
from the request of constant electrochemical potential, 
and the relation ( [Ti] ) is used to evaluate the charge trans- 
fer. 
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III. EQUATIONS OF MOTION 

The total free energy is the sum of the components 
discussed above and the Helmholtz free energy of the 
magnetic field, Tm — J drdz^(B — B ) 2 . Note that 
the magnetic energy does not have a layered structure 
but it is given by a 3D integration. The magnetic field is 
the 3D rotation of the vector potential, B = rotA. 

The total free energy, T = J-kin + J~ con + + , is a 
functional of the wave function tp, the vector potential A 
and the normal density n n or- By variations with respect 
to these functions one arrives at equations of motion. 
For the conventional superconductor, this procedure is 
described in Ref. [l2|. We do not present its straightfor- 
ward modification for the layered structure but discuss 
directly the results. 



A. Maxwell equation 

The functional variation 5T /5A yields the Maxwell 
equation, 



V 2 A 



(15) 



The right hand side includes the in-plane 2D currents, 



j$ = — Re$ (-iftV - e*A^-) ty. 



(16) 



In principle, Josephson currents in the c direction may 
be also present. From the translational and mirror sym- 
metries it follows that no Josephson currents flow be- 
tween two neighboring CuC>2 planes. Assuming very 
small Josephson coupling between the chains and planes 
we neglect also the current between these layers. 
The Maxwell equation (15) has the integral form 



i 2 + «, - z y 



(17) 



For z = Zj the vector potential in the Fourier represen- 
tation simplifies to 



A.(k) = ^E^'( k )- 



D 



(18) 



For the YBCO with its large London penetration 
depth, Aloii = 14000 A, the vector potential changes only 
negligibly along the c axis. Indeed, the applied magnetic 
field 9.4 T results in a vortex distance of 157 A, which 
is small compared to Alou ■ The space modulation of the 
magnetic field is thus covered by the lowest Fourier com- 
ponents of the Abrikosov lattice. Since the vortex dis- 
tance is large on the scale of the YBCO lattice constant 
D, one can take the long wave length limit, k — > 0, in 
which ( |lq ) reads 



A(k) = gj(k), j(k) = ~5>(k). (19) 



In the approximation ( |l9|) the vector potential is identical 
in chain and plane layers so that we can skip its index. 



B. Schrodinger equation 

The Schrodinger equation of a system consisting of 
three periodically repeated layers has a complicated no- 
tation. We first write down this equation in its general 
form and then reduce it by symmetries and the approx- 
imation of weak Josephson coupling between chains and 
planes. The result oriented reader can skip reading de- 
tails of the equations in this section. 

The variation 8T /d^ results in the set of Lawrence- 
Doniach equations 



2m 



0. 



j pp (4 'j - 4j) + Jc P - ^--i) 

(-»RV-e*A) 2 ^-+x7^-=0, 



2m 

Jc P - ^7 +1 ) + Jc P - ^ + ) 

+ -^HW-e*A) 2 ^ + X ^ -0. (20) 
2m* J J J 

We have used zero phases of the Josephson terms. 

For B || c, all Cu02 planes are identical, therefore, 
%l>7 = t/jj = i'J+i = i'p- The chain layers are also identi- 
cal, — ipj-i = V'c- This simplifies the set ( ^p|) to two 
2D Schrodinger equations 

1 2 

J cp O/V - ipc) + tt^ (-«^V - e*A) tp p + x P 4>p = 0, 

2J cp (ipe ~ %) + (~ lhV ~ e * A ) 2 + Xc4>c = 0. 

(21) 

The effective potentials, 



_ 2 £con IpT 



2n c 



1 - 



2|</> £ 



(22) 



follow from the variation of the condensation energy. The 
potential in chains consists of the entropy term only. 

If one neglects the effect of the charge transfer on the 
material parameters (n p , n c , e con , J p , 7c), the set of equa- 
tions ([l9]) and ( pl[ - |2^ ) is closed. As in the ordinary GL 
theory this set describes the magnetic properties of the 
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system. In fact, the transferee! density of holes is very 
small compared to the total densities n p and n c . There- 
fore in the first step, the magnetic field B and the wave 
function ip can be evaluated from this set with no regards 
to the electrostatic phenomena. The charge transfer is 
evaluated in the second step from the wave function ip. 



C. Electrostatic potential 

From variations of the free energy with respect to 
the normal densities n^ T one finds the electrostatic 
potential^ 



\ip P 



etp c 



IV>C 



2 dn v 
2 dn c 



2|Vvl 2 



de coll 2\ib p 
dn v n v 



2|^ c 



(23) 



The right hand sides represent local changes of the chem- 
ical potential. Equations (|2^) thus express that the elec- 
trochemical potential remains constant. 

One can see that the electrostatic potential in planes 
differs from the potential in chains. At first glance, these 
potentials differ by the term 



dn v 4 drir, 2 dn v 



(24) 



We have used the Gorter-Casimir formula for the conden- 
sation energy, e con = jj p T^ and the assumption that the 
thermodynamic properties of the superconducting phase 
are dominated by planes. 

For YBCO the term is much larger than any other 
term in the set (p3|). This is due to a rather strong de- 
pendence of the critical temperature T c op, the density 
of holes in planes. Blatter's approximationEJ of the elec- 
trostatic potential includes the second term of (|^) only. 
For simplicity we use Blatter's approximation by taking 



7 P T C dT c \ip p \ 
e dn„ n n 



if c = 0. 



(25) 



Using approximation (|2^) in the Poisson equation ( |l4| ) 
one obtains the 2D density of holes needed to evaluate 
the electric field gradient effect on the NMR lines. 



D. Limit of weak Josephson coupling 

Within the approximation (^) the electrostatic poten- 
tial depends only on the wave function in planes. For a 
weak coupling between chains and planes, J cp — > 0, the 
set of two non-linear equations (|2l]) simplifies to a single 
equation, 



1 



2m? 



(-iW 



*A) 2 v + x P ^ = o. 



(26) 



We have introduced a 3D wave function if> 



D 'rPi 



terms of which the current j has the usual GL form 

j = — ReiM-i7iV-e*A)V>. (27) 
m* 

The wave function ip corresponds to the 3D density of 
pairable holes, n — 2n p /D. As one expects, in the limit 
of weak Josephson coupling, the magnetic properties are 
described by the customary GL theory. 

To summarize our approach, we use the numericaLcpde 
developed originally as a solver of the GL equations .lia Its 
modification to the Bardeen's approximation has been 
discussed in Ref. |l2|. The output is the magnetic field 
B and the wave function ip. The wave function is scaled 
to its layered counterpart ip p and used in ( p5| ) to provide 
us with the electrostatic potential in planes. From the 
Poisson equation (|l4j) we obtain the density of holes in 
the Cu0 2 planes as 



2fce(l 



-kD\ 



( 1 _ e -2fcD c )( 1 + c -feA»t) 



<p p (k). 



(28) 



In our numerical treatment below we use the full ex- 
pression (p8|). It should be noted, however, that the long 
wave length limit k — > 0, i.e., p P (k) = -p-<^ p (k), gives 
nearly the same result. As the relative errors due to the 
long wave length approximation are of the order of 10~ 3 , 
one can say that the charge is transferred locally from 
planes to chains. If the charge of planes and chains in 
the vortex core is summed together, the vortex cores re- 
main nearly neutral. 



IV. ELECTRIC FIELD GRADIENT EFFECT 

Now we are ready to evaluate the NMR lines. The 
energy levels of the nucleus with spin / = 3/2 
and spin component along the magnetic field, m — 
-3/2,-1/2,1/2,3/2, readES 



E m = -h'fcuBm ■ 



e 2 gQ 
12 



1 3 • 2 

1 sin 

2 



o 2 15 

3m 

4 



(29) 



Here Q is the quadrupole moment of the nucleus, q is the 
EFG which depends on the transfered charge, and 9 is 
the angle between the principal axis of the EFG tensor 
and the magnetic field. 

The NMR frequencies are differences between neigh- 
boring levels, Uvi = E_ X ji - E 1/2 , hu 2 — E 1/2 - E 3/2 , 
and hv-i = E_ 3 /2 — E-i/%. The empirical formula (|l|) 
corresponds to the magnetic field parallel to the axis c. 
When the magnetic field declines from the axis c by angle 
0, the EFG effect on the NMR frequencies is reduced, 



^2,3 = ICuB t (C + AN) 1 



^0 



(30) 
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The amplitude of the local magnetic field changes in 
space, B = (B) + AB. Although the deviations AB from 
the mean field are rather small, they bring an apprecia- 
ble contribution to the NMR frequencies on the scale of 
the shift Az/. The EFG effect depends on the angle 9 
given by disorientations of grains in the sample. In the 
appendix we show that the contribution of the magnetic 
field perpendicular to the axis c which appears even in 
a perfectly oriented crystal due to the layered nature of 
diamagnetic currents is negligible. 

For a quantitative discussion of the EFG effect, it is 
advantageous to compare the NMR frequency with its 
value in the absence of diamagnetic currents, 2/2,3 — i4 3+ 
Av% 3, where 



4% = lCu(B) T AN 



The deviation caused by the diamagnetic currents and 
the disorientation is 



Ai/ 2 , 3 = ICuAB t (aAN - i/q I sin 2 1 



The density of holes per Cu site needed in (|35 
portional to the charge density in the planes 



AN(r) = 



n 



Cu 



(32) 



is pro- 



(33) 



where ficu = 14.88 A 2 is the area per Cu atom in a single 
plane. The density p® describes the transfer of holes in 
the non-magnetic state. Since the non-magnetic state is 
isotropic, i.e., it has only the component k — 0, from (I2F 
and (By) one finds 



Pi 



e \i>o 



D 



-l P Tl 



,<91nT c 
d In n D 



(34) 



where iV'ool 2 = 5(1 — t 4 )n p is the non-magnetic value of 
the GL function in Bardeen's model, t = T/T c . 



f(9) oc sin(#) exp 



k B T, 



B prep 



(36) 



For E OT ^> /csTprop, which is a necessary condition to 
prepare a well oriented sample, the large angles are very 
unlikely and one can use the approximation 



f(0) = 7r 9e- 



2 /e 



(37) 



where Q 2 C = kBT prep /E or measures the scatter of angles. 

The average over the orientation of grains has to be 
performed in addition to the natural line width. The 
observed distribution of frequencies thus reads, 



(31) F 2>3 (u) = / " d0f(0) / dr 



{v- 1/2,3(0, r)) 2 +T 2 ' 



(38) 



where z/2,3 is given by m and @. 

The integral over 9 in d38| ) can be evaluated analyti- 
cally in terms of the exponential integral of complex argu- 
ments. This step, however, is not favorable numerically. 
A faster numerical scheme is obtained if one first evalu- 
ate the distribution F 2; 3 of an ideally oriented sample via 
(||) with the NMR frequency of a perfectly aligned mag- 
netic field (0). The angular averaging is then covered by 
a convolution, 



dxe~ x F. 



2,3 



(39) 



From the experimental data of KNM we estimate T ss 
140 kHz and §z/ 6> 2 = 0.9r = 126 kHz. For these values 
the maximum of the lower NMR frequency is shifted by 
78.3 kHz while the maximum of the upper frequency is 
shifted by —78.3 kHz. This shift appears at any temper- 
ature and it is subtracted in the evaluation of the EFG 
effect. The upper integration limit, 7r 2 /46' 2 s» 10 3 , can be 
replaced by infinity or a smaller cut off at convenience. 



A. Average over grain orientation 

So far we have discussed a single crystal. The sample 
measured by KNMEj is made, however, from c oriented 
grains. Of course, there is a small but finite scatter in the 
orientation of individual grains. The angle 9 then varies 
from grain to grain keeping its value independent of the 
temperature and the magnetic field, provided that the 
applied magnetic field is parallel to the mean orientation. 
The observed NMR lines are the average over all grains. 

We assume that the azimuthal and the polar angles 9 
and 4> are given by the Boltzmann distribution, 



f(9, <j>) cx exp - 



E n 



k B T, 



■ sm 



prep 



(35) 



where E m is the energy needed to disorient the grain 
during a preparation at temperature T prep . One can in- 
tegrate out the polar angle, 



V. NUMERICAL RESULTS 

Our numerical treatment follows the three steps men- 
tioned above. In the first step we evaluate the GL func- 
tion ^i(r) and the magnetic field B(r) of the Abrikosov 
vortex lattice using Bardeen's extension of the GL the- 
ory represented by the set of equations (|10,p*27). In 



this step we take the magnetic field as parallel to the 
axis c neglecting its small perpendicular component in 
individual grains. 

In the second step we use Blatter 's approximation ( |25| ) 
to evaluate the electrostatic potential in planes from the 
wave function. The charge density is obtained from the 
Poisson equation (|2|) and is rescaled to the transfered 
number N(r) of holes per Cu site via (33). 



In the third step we evaluate the local NMR frequen- 
cies for the magnetic field parallel to the axis c using (|3j3) 
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with 9 = 0. The distribution F of the NMR frequencies 
for a single crystal with B || c is obtained via space av- 
eraging of Lorentzian lines (||). The distribution F of 
the NMR frequencies for the granular sample is obtained 
from F via the convolution (j39|). The maximum frequen- 
cies ^2,3 of the crystal and granular samples are found 
from the maxima of -£2,3 and i*2,3, respectively. The lat- 
ter one can be compared with the experimental results 
of KNM, the former represents our prediction for simi- 
lar measurements that, as we hope, will be performed on 
single crystals in future. 

The material parameters of the YBa2Cu307 from 
Rcf. |2l] are: the density of pairable holes n = 5 x 
10 27 m -3 , the effective mass m* = 6.92 m e , the critical 
temperature T c = 90 K, and the Sommerfeld constant 
7 P = 302 JK -2 m -3 . These values give a London pen- 
etration depth Alod = 1-4 x 10~ 7 m, an upper critical 
field B C 2 = 96.5 T, and a GL parameter at the criti- 
cal temperature -kq = 55. In Fig. 2.16 PlakidaEj shows 
results of JunodE 2 ] according to which the charge trans- 
fer — 0.03e from chains to planes per Cu site leads to a 
decrease of the critical temperature by 30 K. This corre- 
sponds to d In T c /d In n p = —4.82. We take the permit- 
tivity e = 4e - We do not discuss the YBa 2 Cu40s for 
which we could not find relevant material parameters. 



A. Space variation of the NMR frequency 

Figure [l] shows the fishnet plot of the lower NMR fre- 
quency At^(r) for T = 0.6 T c and 9 = 0. The variation 
is of the order of 100 kHz, which is comparable to the 
line width. 
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FIG. 1. The lower NMR frequency Ai^(r) in the crystal at 
T = 0.6 T c . The frequency reflects the triangular structure of 
the Abrikosov lattice. The maxima are at the vortex cores. 

In Fig. ||we plot the contribution of the Zeeman effect, 
Ai^z = iCuAB, and the EFG effect, Ai/ E fg = —AAN 
separately. One can see that the Zeeman effect is about 
seven times larger that the EFG effect. The main con- 
tribution to the spatial variation of the NMR frequency 



is thus due to the compression of the magnetic field in 
vortex cores. On the other hand, the EFG has nonzero 
mean values, therefore, it contributes more to the observ- 
able shift of the maxima of the NMR line. 
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FIG. 2. The Zeeman effect Au z (r) = 7c„AB(r) (upper), 
and the EFG effect Ai/efgM = -AAjV(r) (lower). Parame- 
ters are the same as in Fig. hi 

At temperatures close to T c , the space variation of the 
Zeeman effect is reduced following the reduced compres- 
sion of the magnetic field in the vortex cores. Fig. || shows 
a decomposition of the lower NMR line into the Zeeman 
part and the EFG effect for T = 0.9 T c . One can see that 
here the amplitude of the EFG effect is close to the am- 
plitude of the Zeeman effect. Note that the mean value 
of the EFG is larger than its spatial variation. This sig- 
nals that the mean value provides a good approximation 
at this temperature. 
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FIG. 3. The Zeeman effect Av z (r) = 7 Cu A_B(r) (upper), 
and the EFG effect Az/efgM = -AAN(r) (lower), for 
T = 0.9 T c . 



B. NMR lines of a single crystal 

Above the critical temperature, the NMR lines of a 
single crystal are given by the Lorentzian distribution. 
Below T c the space variation of the magnetic field de- 
forms the lines so that their centers of mass (i.e., mean 
values) are not identical to their maxima, see Fig. 4. 

The deformation of the NMR line seen in Fig. 4 can 
be correlated with the space distribution of the NMR 
frequency presented in Fig. [I]. The vortex cores supply 
frequencies in the range from zero to 75 kHz above the 
mean value. Their contribution is visible as the extended 
left shoulder of the line in Fig. ||. The intermediate region 
between vortices supplies frequencies about 25 kHz below 
the mean value. Their contribution results in a shift of 
the maximum of the NMR line. 



v [kHz] 

FIG. 4. The distribution F2(v) (solid line) compared with 
the Lorentz distribution F^, OT {y) — ^ ly a^ r a (dashed line) at 
T = 0.6 T c and Y = 50 kHz. The experimental line width is 
F = 140 kHz. We use the narrow line to make its deformation 
more visible. 

The temperature dependence of the mean values 
(Ai/2.3) and the maxima Ai/2,3 is presented in Fig. |5[ 
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FIG. 5. The maxima £"2,3 (dashed lines) and mean values 
(A^2,3) (full lines) of the lower and upper NMR lines (be- 
low) as a function of temperature. The corresponding shifts, 
Av = \ (Au 2 ~ Av 3 ) and (Av) = |((Ai/ 3 ) - {Au 2 )) are plot- 
ted above. 

The mean value is independent of the space variation 
of the magnetic field, (Az/2,3) = 7c u (A J B) =F A(AN) = 
^A(AN), because (AB) = by definition. The differ- 
ence between the mean value and maximum thus shows 
the effect of the spatial variation of the magnetic field 
on the position of the maximum. Since the mean value 
is proportional to the charge transfer while the position 
of the maximum is a rather complicated quantity, the 
temperature region in the vicinity of T c is more recom- 
mendable for the experimental exploration. We reserve 
the symbol T c for the critical temperature in the absence 
of the magnetic field, T c = 90 K. In the magnetic field 
B = 9.4 T, the actual phase transition is at T = 0.931 T c . 



The deviation of the maximum from the mean value 
depends on the line width. This is illustrated in Fig. ^ 
by the line of artificially chosen width T = 50 kHz. In 
the limit of very broad lines, r — > oo, the position of 
the maximum becomes identical to the mean value. This 
trend is also seen from the formula (0). 



C. Conventional model 

Let us make a small detour and assume that there is 
no charge transfer between chains and planes. In other 
words, no screening by chains is assumed and the charge 
is transfered exclusively inside the CuC>2 planes. In this 
case the system would behave as a conventional supercon- 
ductor in which the charge is transfered from the vortex 
cores into the region between them. In this model, the 
charge density is given by p p = ip p /L pp , while other rela- 
tion remain the same as in the case of the above layered 
model. 

One can see in Fig. || that the sign of Av is oppo- 
site to the sign of the model with screening by chains. 
Also, the magnitude of Aj/ is extremely small. A typi- 
cal shift is less than 1 Hz while the typical shift is more 
then 1 kHz in the case where the chains are accounted 
for. The values observed by KNM are about 10 kHz. 
Our result thus confirms the estimate of KNM that the 
conventional model of the. vortex charge fails to explain 
their experimental data.tU 

When the transfer between chains and planes is prohib- 
ited, the amplitude is reduced by two mechanisms. First, 
the charge transfer accross the vortex core happens on a 
distance which is much larger than the distance between 
chains and planes. Second, the mean value of the shift is 
zero, (i>i) = 0, because the charge of planes is conserved, 
Nq — (N) = 0. The line shift thus results merely from 
the statistical variations which are small due to the large 
line width, see (|J). 

An explanation of the reversed sign follows the argu- 
ment given already by KNM.E2I Since the magnetic field is 
much smaller than H C 2 , the area of vortex cores is much 
smaller than the area between them. Accordingly, cores 
have a negligible weight in statistical variations so that 
the charge between cores determines the line shift. This 
charge is opposite to the charge of cores, consequently 
the reversed sign of the line shift results from the con- 
ventional model. 
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FIG. 6. The maxima i>2,3 of the lower and upper NMR lines 
as a function of temperature (below) and the shift Av (above) 
in the absence of screening chains. 



D. Averaging over grain orientation 

As mentioned, KNM have used a granular sample, with 
a high but not perfect orientation of YBCO grains .E3 The 
averaging over the grain orientation ( j39|) results in an ad- 
ditional shift of the NMR frequencies, see Fig. @. The 
grain disorientation contributes mainly at low tempera- 
tures where it leads to shifts of about 5 kHz, while the 
value predicted for the single crystal reduces with tem- 
perature. 
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FIG. 7. The maxima Ai/2,3 of the lower and upper NMR 
lines as a function of temperature (below) and the shift (Au) 
(above) of T = 140 kHz of Fig. |E] compared with the averaging 
over grain boundaries according to (B9h . 



Figure g compares our theoretical results with the ex- 
perimental data of KNM. Apparently, the agreement is 







only qualitative. We have obtained the correct sign but 
the theoretical amplitude is about two or three times 
smaller than observed. 
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FIG. 8. Comparison of the theory with experiment. The 
experimental data are the dots with error bars, the theoretical 
prediction is the full line and the theoretical value for the 
single crystal is shown as a dashed line. 

In the temperature dependence one can identify three 
regions. Above the critical temperature there is no charge 
transfer between chains and planes. Indeed, neglecting 



the second terms in the potential (23) we have neglected 



the thermally induced charge transfer in the normal state 
which appears for unequal energy derivatives of the den- 
sity of states_97 P /9n p ^ dj c /dn c . The experimental 
data of KNMEj above the critical temperature justify this 
neglect. 

In the narrow region between T c and 0.93 T c the sam- 
ple is superconducting only if the magnetic field is absent. 
The NQR frequency is thus measured in the supercon- 
ducting state while the NMR frequencies arc measured 
in the normal state. In this temperature region the ob- 
served frequencies are free of space fluctuations giving a 
direct experimental access to the charge transfer caused 
by the superconducting transition. Since in both cases 
(NMR and NQR) the magnetic field and the transfered 
charge are homogeneous, the averaging over the grain 
orientation is the same as in the normal state above T c . 
Accordingly, the granular sample used by KNM should 
provide data comparing well with the single crystal. Un- 
fortunately, this region is not covered by sufficiently de- 
tailed data in Ref. |l3|. 

The region between 0.93 T c and T c can be used to ex- 
perimentally test our assumption about the charge trans- 
fer between chains and planes. In the absence of the 
charge transfer between chains and planes, the shift Av 
is exactly zero since there are no vortices to cause the 
in-plane transfer. 

Below 0.93 T c the sample is superconducting even in 
the presence of the magnetic field, i.e., the Abrikosov 
vortex lattice enters the game. Let us first discuss the 



mean value. When the magnetic field penetrates the 
sample, the superconductivity is suppressed in the vor- 
tex cores. The charge transfer induced by vortices is 
crudely proportional to the product of the superconduct- 
ing fraction and the area of the vortex core. As the tem- 
perature is reduced, the superconducting fraction grows 
while the diameter of the vortex cores shrinks down. For 
a GL parameter re independent of the temperature, the 
product remains constant. Within Bardeen's approxi- 
mation the GL parameter re slightly decreases with the 
temperaturej!3 therefore, the charge transfer reduces as 
T ^> 0. 

At temperatures below 0.7 T c the maximum of the 
NMR line in the granular sample is strongly shifted by 
the angular averaging. According to our results, the 
granular samples in this temperature region cannot be 
used to access material properties. 



VI. SUMMARY 

We have shown that the experimental results of Ku- 
magai, Nozaki and Matsuda can be explained in terms of 
a standard theory of the electrostatic field in supercon- 
ductors extended by features necessary to describe the 
layered structure of YBCO. The agreement of our theory 
with the experiment is only qualitative. Our results, thus 
do not exclude possible alternative mechanismsEJ'EJ. 

In conclusion, the charge transfer between CuO chains 
and Cu02 planes in YBa2Cu307 provides a qualitative 
explanation of the electric field gradient effect on the 
NMR lines. At low temperatures the effect of the charge 
transfer is spoiled by two additional effects: the space 
modulation of the Zeeman effect reduces the shift of the 
maxima of the NMR lines; the disorientation of grains in 
the sample increases the shift. From this point of view it 
is advisable to make precise measurement near the criti- 
cal temperature where the charge transfer dominates. 

Our study does not cover two realistic features that 
can increase the amplitude of the observed shift of NMR 
lines towards the experimental values. First, the charge 
is transfered from the copper to the apex oxygen, i.e., the 
distance between planes and chains D cp = 4.24 A should 
be reduced to the Cu-0 distance 2.3 A. This will enhance 
the charge transfer by a factor 1.8. We did not include 
this simple correction as its consistent implementation 
requires to use the non-local dielectric function and the 
microscopic treatment of electronic states. A similar type 
of enhancement would result from a higher value of the 
dielectric constant. We have used e = 4e r*jhile val- 
ues as large as e = 100eo has been advocated.c3 Second, 
the magnetic anisotropy of YBCO enhances the devia- 
tions of the magnetic field from the axes c in individual 
grains. At low temperatures this effect should be stronger 
in agreement with the observed trends. This correction 
is omitted because of an unknown vortex pinning. 

It is apparent from the experimental data of KNM 
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that an explanation of the NMR lines in the underdoped 
YBa2Cu408 will be more complicated. Unfortunately, 
we could not discuss YBa 2 Cu 4 08, because we have not 
found a relevant set of material parameters. 
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APPENDIX A: INNATE PERPENDICULAR 
MAGNETIC FIELD 

The magnetic field perpendicular to the axis c appears 
in the sample from two sources. First, the observed sam- 
ples are made from many grains with a small but finite 
scatter of the orientation of the c axis. This contribution 
is discussed in the main text. Second, since diamagnetic 
currents are in planes, the magnetic field is more com- 
pressed in layers of planes than in layers of chains. This 
part of the perpendicular field is innate to the layered 
system and appears also in single crystals. 

Here we show that the innate perpendicular magnetic 
field can be neglected. The magnetic field is given by the 
vector potential A. For B || c, the component of A along 
the c axis is zero, see (|i~7j). The perpendicular magnetic 
field thus reads, 



B , 



dAy 

dz 



dA fl 



By — 

O 



(Al) 



z x <9 Z A, where z is the unit 



In the short notation, Bj 
vector in the direction z. 

The vector potential ([l7]) has a cusp at the positions of 
layers. This cusp appears due to the singular 2D current 
of the selected layer, say at zj~ . From the mirror symme- 
try follows that the current in the z^-layer does not con- 
tribute to the perpendicular magnetic field in the center 
of this layer, i.e., at z+ where the Cu nuclei sit. We can 
thus avoid the cusp by taking neighbor layers only. The 
mirror symmetry also shows that all equivalent z~p layers 
spaced symmetrically above and below z~^ give zero net 
contribution. Accordingly, the perpendicular magnetic 
field in the vicinity of Zq (for \z — D pp /2\ <C D pp /2) is 
given by the reduced vector potential, 



A Iea (r,z) = — I dr 
4ir 



|r' - r| 2 + (zj, - z) 2 



(A2) 



The 2D Fourier transformation of (A2) yields 



A red (k,z) = 



Mo e 



-k(z+D pp /2) 



-k{D-z-D pp /2) 



2k 



1 



-kD 



-jp(k). 

(A3) 



The derivative of A with respect to z at z — D pp /2 
provides us with the perpendicular magnetic field in 
terms of the 3D current, j = 2] p / D, 



B 



HoD e 



kD„ 



-k(D-D pp ) 



1 



-kD 



-zxj(k). (A4) 



Let us estimate the value of the perpendicular field. To 
this end we employ the long wave length limit, k — > 0, 



B 



Mo 



(D - 2D w )i x j(k). 



(A5) 



The maximum of the current is at the edge of the vortex 
core, i.e., at r w v2£j where £ = X/k is the GL coherence 
length. The amplitude of the vector potential at this 
point, A = $o/(27rr) = &ok/(2it\/2\), gives the current 
density fj, j = A/X 2 = <P k/(2^ttX 3 ). For the YBCO 
parameters, [D — 2D pp )/4 — 1.33 A, k = 55, and A = 
1400 A, one finds B ± ^6x 10~ 4 T. This field is far too 
small shifting the NMR line by i/ | sin 2 6 w 0.2 Hz. We 
have used v = 31.5 MHz. 



1 G. Rickayzen, J. Phys. C 2, 1334 (1969). 

2 J. Bok and J. Klein, Phys. Rev. Lett. 20, 660 (1968). 

3 J. B. Brown and T. D. Morris, Proc. 11th Int. Conf. Low. 
Temp. Phys., Vol. 2, 768 (St. Andrews, 1968). 

4 T. D. Morris and J. B. Brown, Physica 55, 760 (1971). 

5 Yu. N. Chiang and O. G. Shevchenko, Fiz. Nizk. Temp. 12, 
816 (1986) [Low. Temp. Phys. 12, 462 (1986)]. 

6 Yu. N. Chiang and O. G. Shevchenko, Fiz. Nizk. Temp. 22, 
669 (1996) Low. Temp. Phys. 22, 513 (1996)]. 

7 P. Lipavsky, J. Kolacek, J. J. Mares and K. Morawetz, 
Phys. Rev. B 65, 012507 (2001). 

8 D. I. Khomskii and F. V. Kusmartsev, Phys. Rev. B 46, 
14245 (1992). 

9 D. I. Khomskii and A. Freimuth, Phys. Rev. Lett. 75, 1384 
(1995). 

G. Blatter, M. Feigel'man, V. Geshkenbein, A. Larkin and 
A. van Otterlo, Phys. Rev. Lett. 77, 566 (1996). 

11 J. Kolacek and P. Lipavsky, Physica C 364, 138 (2001). 

12 P. Lipavsky, J. Kolacek, K. Morawetz and E. H. Brandt, 
Phys. Rev. B 65, 144511 (2002). 

13 K. I. Kumagai, K. Nozaki and Y. Matsuda, Phys. Rev. B 
63, 144502 (2001). 

14 H. Yasuoka, in Spectroscopy of Mott Insulator and Corre- 
lated Metals, editted by A. Fujimori and Y. Tokura, Solis 
State Sciences 119, 213 (Springer- Verlag, Berlin 1995). 

15 W. E. Lawrence and S. Doniach, in Proc. of the 12th Int. 
Conf. on Low- Temperature Physics, Kyoto, 1970, edited by 
E. Kanda (Keigaku, Tokyo, 1971), p. 361. 



11 



I. N. Bulaevskii, Int. J. of Modern Physics B 4, 1849 (1990). 
J. Bardeon, Phys. Rev. 94, 554 (1954). 
J. Bardeen, Theory of Surperconductivity in Handbuch der 
Physik, Bd. XV. (1955). 

E. H. Brandt, Phys. Rev. Lett. 78, 2208 (1997). 

C. P. Slichter, Principles of Magnetic Resonance (Springer- 

Verlag, Berlin 1978). 

N. M. Plakida, High- Temperature Superconductivity 
( Springer- Verlag, Berlin 1995). 

A. Junod, in Physical Properties of High Temperature Su- 
perconductors II, ed. D. M. Ginsberg (World Scientific, Sin- 
gapore 1990). 

E. Simanek, Phys. Rev. B 65, 184524 (2002). 
Y. Chen, Z. D. Wang, J.-X. Wang and C. S. Ting, The 2002 
Int. Conf. on Phys. and Chem. of Molecular and Oxide Su- 
perconductors (MOS2002), August 13-18, 2002 in Hsinchu, 
Taiwan. To appear in J. Low. Temp. Phys. 
J. Mannhart, Mod. Phys. Lett. B 6, 555 (1992). 



12 



